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We investigate a semiconductor p-n junction in contact with superconducting leads that is oper¬ 
ated under forward bias as a light-emitting diode. The presence of superconductivity results in a 
significant increase of the electroluminescence in a sharp frequency window. We demonstrate that 
the tunneling of Cooper pairs induces an additional luminescence peak on resonance. There is a 
transfer of superconducting to photonic coherence that results in the emission of entangled photon 
pairs and squeezing of the fluctuations in the quadrature amplitudes of the emitted light. We show 
that the squeezing angle can be electrically manipulated by changing the relative phase of the order 
parameters in the superconductors. We finally derive the conditions for lasing in the system and 
show that the laser threshold is reduced due to superconductivity. This reveals how the macroscopic 
coherence of a superconductor can be used to control the properties of light. 
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I. INTRODUCTION 

Superconductors exhibit quantum coherence of elec¬ 
tronic degrees of freedom on a macroscopic scale. For 
s-wave pairing they are characterized by a complex pair¬ 
ing amplitude A with a well-defined phase . This is a 
consequence of the off-diagonal long-range order present 
in the two-particle density matrix^’^. The electrons con¬ 
dense into Cooper pairs and form entangled two-electron 
singlet states. The prospect of harvesting these useful 
coherence and entanglement properties for the manip¬ 
ulation of quantum states as required for example in 
quantum information processing and communication^^® 
is a motivation to integrate superconducting elements 
in semiconductor solid-state nanostructures. Prominent 
specific goals are the on-demand production of entan¬ 
gled photon pairs due to the recombination of Cooper 
pairs^^^® and the generation of non-classical states of 
light^^“^®. It further brings about the fundamental ques¬ 
tion of how to efficiently transfer the electronic coherence 
and entanglement naturally present in a superconductor 
to excitonic particle-hole pairs in a semiconductor and 
eventually to photons that are emitted from the hetero¬ 
structure. 

While the experimental realization of superconductor- 
semiconductor hybrid nanostructures has proven to be 
technology challenging^^, there have recently been a 
number of successful experiments using semiconduc¬ 
tor nanowires^®’^®, quantum wells®’®® and self-organized 
hetero-superstructures in unconventional superconduc¬ 
tors®^. Various desired electronic and optoelectronic 
properties have been observed such as the proximity ef¬ 
fect®’®® and the Josephson effect^®’®®, the realization of 
a superconducting field-effect transistor®® and enhanced 
electro- and photo-luminescence rates®’®®’®^^®®’. 

Enhanced luminescence rates were observed in a se¬ 
ries of experiments using superconductor-semiconductor 
hetero-structures where a p-n junction was contacted 
with a superconducting Niobium lead on the n-side of 
the junction and cooled below the superconducting tran¬ 


sition temperature X),®’®®’®^’®®. In addition, it was clearly 
shown via the Josephson effect that due to the prox¬ 
imity effect Cooper pairs tunnel into the active region 
of the p-n junction, where they recombine with normal 
holes®. More recently, enhanced photo-luminescence be- 



FIG. 1. (color online). Schematic setup of superconduct¬ 
ing light emitting diode with p-n junction coupled to su¬ 
perconducting (sc) leads operated under forward bias volt¬ 
age Vq. Basic recombination processes at low temperature 
T < |Ac|,|A„|: (i) transfer of particle from conduction to 
valence band upon radiation of a photon (red) with energy 
ujq < eVb — I Ac I — |A„|, (ii) absorption of a photon (blue) 
with energy uoq > eVb -|- |Ac| -I- |Ac| upon transfer of a par¬ 
ticle from valence to conduction band, and (iii) Cooper pair 
tunneling upon emission or absorption of two photons (green) 
with energy ujq — eVb- The superconducting gaps are denoted 
Ai, (Ac) for valence (conduction) band. 
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low a superconducting transition temperature was re¬ 
ported in self-organized superconductor-semiconductor 
hetero-structures that form naturallly in the iron-based 
superconductor Ka;Fe 2 -ySe 2 ^^. The exciting case where 
the semiconducting region only contains a few quantum 
dot states has been studied in Refs. 10, 11, 28-30. 

Here, we investigate a superconductor-semiconductor 
hetero-structure consisting of a p-n junction that is sand¬ 
wiched between two superconducting leads. The setup 
is shown in Fig. 1. We focus on the properties of the 
light that is emitted from such a superconducting light- 
emitting diode (LED) under forward bias in the steady- 
state. To reach a steady-state we consider a coupling of 
the photons to an external photon bath. One of our 
main results is a drastic enhancement of the electro¬ 
luminescence rate in a sharp spectral window in the pres¬ 
ence of superconductivity. This increase, which is shown 
in Figs. 6 and 7 is a consequence of the van-Hove singular¬ 
ities at the edges of the superconducting bands. Similar 
behavior was observed in a setup with one normal and 
one superconducting lead in Refs. 26 and 27. 

We also report the formation of an additional lumines¬ 
cence peak on resonance that is effectively due to transi¬ 
tions of Cooper pairs from the conduction to the valence 
band. A similar rich peak structure was reported for the 
case of quantum dots coupled to superconducting leads 
in Ref. 10. Further, we analyze the statistical properties 
of the emitted photons and predict non-classical, two¬ 
mode squeezed states of light. Photons with momentum 
q and —q inherit their squeezing correlations from the 
electrons forming Cooper pairs. As shown in Fig. 10, the 
uncertainty of the corresponding two-mode quadrature 
operators in the squeezed state falls below the minimum 
uncertainty of coherent states enforced by the Heisen¬ 
berg limit. Even more importantly, the orientation of 
the squeezing ellipse can be externally manipulated by 
changing the relative phase between the two supercon¬ 
ductors that are coupled to the p-n junction^®. This 
shows that one may employ the coherence of a supercon¬ 
ductor to tailor the properties of a two-photon quantum 
state via the transfer of Cooper pair entanglement to the 
photons. Two-photon correlations in related models are 
theoretically discussed in Refs. 10, 12, and 13. Einally, 
we discuss conditions for lasing in this system and find 
that the lasing threshold is reduced in the presence of su¬ 
perconductivity, again due to the formation of van-Hove 
singularities. 

The structure of the remainder of the paper is 
as follows: in the next Sec. H, we introduce the 
superconductor-p-n-superconductor setup and develop 
the many-body non-equilibrium field theory framework 
that we use to extract properties of the photons that are 
emitted in such a device under forward bias. In Sec. HI, 
we derive the effective photonic action by integrating over 
the electronic degrees of freedom. We calculate the pho¬ 
tonic self-energy due to the coupling to the electrons up 
to one loop. We discuss in detail the structure and the 
features of the self-energy both in the presence and in the 


absence of superconductivity. To obtain a steady state, 
we consider the leakage of photons out of the system via 
a coupling to an external photon bath. This yields an¬ 
other contribution to the photonic self-energy. Inverting 
the Dyson equation, we sum up the complete RPA series 
and obtain our main analytical result: the dressed pho¬ 
ton Keldysh propagators. In the following sections, we 
extract the physical information they contain about the 
photonic system. In Sec. IV, we discuss the luminescence 
and the light squeezing properties and in Sec. V we de¬ 
rive the conditions to obtain lasing in the system. We 
determine the lasing threshold also in a physically very 
transparent way using rate equations. To focus on our 
main results and the physical consequences in the main 
text, we shift details of the calculations that yield the 
effective photon action, the photonic self-energy and the 
dressed photon propagators to the Appendix. 

II. HAMILTONIAN AND KELDYSH ACTION 
OF SUPERCONDUCTING LED 

A. Hamiltonian of the hetero-structure 

We consider a p-n junction that is coupled on each side 
to a superconducting lead. The system is set under an 
external forward bias voltage Vq to operate as a LED. 
Due to the proximity effect Cooper pairs tunnel into the 
active region of the junction, and participate in recom¬ 
bination processes that lead to the emission of light. We 
model the active region of the junction together with the 
photonic degrees of freedom by the Hamiltonian 

H = H, + H^ + Hei-ph + Hph + ■ ( 1 ) 

The Hamiltonian of the conduction and the valence elec¬ 
tron bands takes the form 

fc,cr=t,4' 

+ -I-h.c.) , (2) 

k 

where creates for a = c{v) an electron with mo¬ 
mentum k and spin a in the conduction (valence) band. 
The band dispersions are given by ea{k) which we will 
later assume to be of the form e^k) = k'^l2nic and 
e„(fc) = k'^l2my with < 0 < nic, be., the conduction 
(valence) band is electron (hole) like. We use h = 1 here 
and in the following. The proximity induced BCS s-wave 
gaps are denoted by and the chemical potentials by 
Pa ■ We have included the chemical potentials into Eq. (2) 
to be able to work with time-independent gap functions 
Aq,. This requires that we measure the electronic ener¬ 
gies with respect to the two different chemical potentials 
Pa in the two respective bands. Their difference is equal 
to the applied bias voltage 

pc Pv — eVo , 


( 3 ) 
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which drives the system out-of-equilibrium. Formally, 
this is achieved by the gauge transformation akcr —> 
cifccr = This changes the dispersions from ea{k) 

fO ^a(^) — ^a(^) l^a- 

In addition, this gauge transformation leads to a time- 
dependence of the electron-photon coupling constant 
9o 9{t) = that appears in the coupling Hamil¬ 

tonian 


Hei-ph = - X! +h-c-) • (4) 

k.k'.cr 

The operator 6^ creates a photon with frequency ojq = 
c\q\. The electron-photon coupling is such that the trans¬ 
fer of each electron from the upper conduction to the 
lower valence band leads to an emission of a single (op¬ 
tical) photon with frequency ujq « eVo. We restrict our¬ 
selves to the case of spin-conserving recombinations. The 
generalization to circularly polarized photons is straight¬ 
forward and does not yield any qualitatively new aspects 
to our model. 

The free photon Hamiltonian is given by 

Hph = '^OJqblbq. (5) 

q 

Since photons are constantly produced in the LED due 
to the presence of the bias voltage, we need to include 
an external photon bath to deposit the energy and reach 
a steady-state. This is described by the bath Hamilto- 
nian^^’^^ 

q q,q' 

The photon bath operator creates a photon with fre¬ 
quency Vq in the bath. The coupling constants between 
system and bath photons are denoted by Xqq'- We will 
later integrate over the bath degrees of freedom leading 
to the possibility to dissipate a (system) photon to the 
external bath and to create a (system) photon from the 
bath. 


B. Action on the closed time contour 


that we are always below the lasing threshold and de¬ 
rive parametric conditions where lasing occurs later (see 
Sec. V). 

In non-equilibrium it is required to begin and end the 
time evolution at the initial state of the system. One 
thus formulates the field theory on a closed time contour 
C such that the action reads^^^^® 



Here, Ck,a,Vk,a denote Grassmann fields and bq,aq de¬ 
note complex fields arising in a path integral formulation 
of the action. One can now clearly observe how the ap¬ 
plied bias voltage appears by performing the gauge trans¬ 
formation introduced above aka —> dfco- = e^^^^aka- The 
time derivative then produces additional terms that con¬ 
tain the chemical potentials such that the energy dis¬ 
persions are measured with respect to the two different 
chemical potentials ea{k) ^a{k) = ea{k) — fjLp and 
9o g{t) = 506 “'^“*- We suppress the tilde notation 
from now on. 

The contour C starts (and ends) at the initial time 
t = 0, where we assume that go = 0 and the system is 
in thermal equilibrium at potentially different temper¬ 
atures for the electronic and photonic subsystems. At 
t > 0 we consider a non-zero electron-photon coupling 
go and a non-zero bias voltage Vq > 0 leading to photon 
production. In the following, we focus on steady-state 
properties of the system at times t, t' ^ ts where G is a 
characteristic time-scale over which transient effects as¬ 
sociated with the switch-on decay. We will always check 
self-consistently that our assumption of a steady-state 
holds. 

It is convenient to introduce fermionic and bosonic 
spinors 



( \ 

V(4,.,)V 



Due to the applied bias voltage the system that we 
are considering is in a non-equilibrium state. In the fol¬ 
lowing we want to set up the required non-equilibrium 
Keldysh field theory formalism that is used to determine 
observable properties of the system such as the electro¬ 
luminescence or the squeezing of the emitted light. The 
small parameter that controls our (infinite-order) pertur¬ 
bative calculation is the electron-photon coupling con¬ 
stant go. More precisely, the small dimensionless param¬ 
eter is given by \9o\'^P f/E p ^ 1 for normal conducting 
leads and |gopgF/|A| ^ I for superconducting leads. 
Here, pp denotes the electronic density of states at the 
leads’ Fermi energy Ep (see Fig. 1). We further assume 


where C = i denotes that the time variable is located 
on the forward (+) or the backward (—) branch of the 
contour C. The A superscript denotes the fermionic com¬ 
bined Nambu-conduction/valence space and the o super¬ 
script describes the photon particle-hole space, in which 
the complex bosons can be described with real spinors, 
see Eq. (20). It is convenient to make a transformation 
from the contour (-1-,—) basis to the RAK (Retarded- 
Advanced-Keldysh) basis. Eor bosons this transforma¬ 
tion to the classical and quantum fields is given by 

¥q‘ = {^+ +^q)/V2 (9) 

4>^ = (<5j-$-)/y2. (10) 
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For fermions we follow Larkin and Ovchinnikov^^ and 
perform the transformation to ( 1 , 2 ) fields as 

( 11 ) 

U = (12) 

A, 

The conjugate Grassmann variable fields (^/fc)' are not 
A , 

related to the 4 /^ fields and one may thus choose a dif¬ 
ferent transformation for them 

(^l)^ = ((fe^ - (^fc )^)/v^ (13) 

(fe^ = ((fe^ + (^fc)^)/V 2 - ( 14 ) 

Combining the two fields as usual we introduce the 

o o » A A 

Keldysh vectors $q = (4)q,<I)^) and 4'^ = (4/fc,4'fe). 
We denote these vectors by A superscripts. We can then 
write the Keldysh action as 

s = r 

J-OO L ^ 

Q 

+E -1')|+^pr ■ (is) 


with dispersions ^aik) = ea{k) — Ha for a = c,v and 
Pauli matrices tTx.y.z where a± = ^(cr^ ± fo^). Since 
the electrons are assumed to be in thermal equilibrium, 
one can express the Keldysh Green’s function via the 
fluctuation-dissipation theorem as 


= F(cv)[G^Jcv)- G^,^H]. (19) 

with fermionic distribution function F(uj) = 1 — 

2np{uj) = tanh(w/2Tp) at temperature Tp. 

The free Green’s function for the photons in the system 
described by bq can be written as 


Do,q{t,t') 


-z($q(t)rq(t'))o = 


0 ) ■ 

( 20 ) 


Note that the expectation value contains the transpose 
spinor 4>!))q(t') (and not the hermitian conjugate). The 
retarded and advanced blocks take the form 


R/A 

0 , 


it-t') 


( 0 

[d^i^it'-t) 0 )■ 


( 21 ) 


where [d^^^( uj)]~^ = uj — ujq ± iO. Since the photonic 
subsystem is initially in equilibrium, we may write the 
Keldysh component via the fluctuation-dissipation theo¬ 
rem as 


The factor of 1/2 in front of the bosonic propagator arises 
from our choice of “real” bosonic spinors. The unper¬ 
turbed fermionic Green’s functions are given by 

GoMbt') = . 

V 0 Go.fefoi')/ 

(16) 


The retarded and advanced blocks, which only depend 
on the time difference t = t — t' , explicitly read 


^R/A 

^0,k 

(t 

-t') = 

■/ ^1/2 


7o 


/ ^{p),R/A 

0 

pR/A 

^0,k,v 

0 



0 

gI^IT 0 

jjRjA 

-^0,fc,c 



^R/A 

^0,k,v 

0 

^{h),R/A 

^0,k,v 

0 


\ 

0 

pRiA 

'^0,A;,c 

0 

Mh),R/A 

^0,k,c 


The average ( • )o is with respect to the free action 
S{go = 0 ). The particle and hole propagators are denoted 
by and The anomalous electronic propaga¬ 

tors which arise due to the presence of superconductivity 
are denoted by P and P. After Fourier transformation 
G(a;) = dTG(T)e^‘^'^ they take the explicit form 


Do^qH = BoH [D^q(uj) - D^q(uj)] . (22) 

Here, Bq denotes an initial bosonic distribution function 
of the uncoupled {go = 0) system. As shown below the 
distribution function Bq will not be important in our 
calculation of the photon distribution in the steady state. 
This is rather determined by the interplay between the 
photon distribution function in the external photon bath 
and the photon production rate in the LED. 

The electron-photon coupling part in the action con¬ 
tains the vertex expression 

14-fe'(t)= E E ^*(^)‘j’^fe'.*(^) (23) 

OL—cl^q i— 1^2 


with matrices 7 “* = 1 and 7 '^ = dx in Keldysh space. 
Using Eq. ( 8 ), the components of the photon field vector 
read explicitly $fc_fc,,i(t) = K-k'i^) and , 2 ^^) = 

(fe“fc+fe/(l))^- The coupling matrices take the form 


A 

9i{t) 


1 

71 


/ 0 0 0 0 \ 

g{t) 0 0 0 1 

0 0 0 -g{t) 

V 0 0 0 0 / 


and 


(24) 


( / \ \ R/A. 

GqI, Po,a \ _ CU ■ 1 + ^a(fc) ■ - Ac ■ g+ - A* ■ g_ 

^o.a ” iu;±^0y-Uk)^-\A^\" 

( 18 ) 


A 

52 (t) 


1 

71 


/o 9{t) 
0 0 
0 0 
Vo 0 


0 0 \ 
0 0 I 

0 0 
-g{t) 0/ 


(25) 
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with g{t) = and g{t) = {g(t))* denoting the com¬ 

plex conjugate. Finally, the part of the action that de¬ 
scribes the coupling of the photons to the photon bath 
reads 


Qbath 

'^ph 


dt 


^a\{idt - Vq)c 


+ 'y + c.c.) 

(1,Q' 


(26) 


We will later in Sec. IIIB integrate over the external 
photon modes {oq} under the standard assumptions of 
an Ohmic bath in the white noise limit of frequency in¬ 
dependent couplings Aqq' and bath density of states^^. 


III. EFFECTIVE PHOTON ACTION 


In this section, we want to derive an effective action for 
the photonic degrees of freedom in the system S'®® which 
takes the coupling to the superconducting leads as well as 
to the photon bath into account. Formally, we integrate 
over the fermion fields and the external bath pho¬ 
ton fields Uq. This integration yields an electronic con¬ 
tribution to the photon self-energy and a bath 

contribution to the self-energy <'). The effective 

photon action takes the form 


CCff _ 

^ph 




POO 


dtdt' (t) 


— OO 


Q 


Do,lit, 


$_q(0. 


(27) 


While we treat the bath self-energy exactly, we take the 
electronic self-energy up to one loop into account. Impor¬ 
tantly, due to the presence of superconductivity there ex¬ 
ist non-zero anomalous fermionic Green’s functions P, P. 
This leads to anomalous terms in the bosonic self-energy 
ri®*(t,t'), which induce similarly anomalous photon ex¬ 
pectation values such as (^q(i)^-q(O)- 

In the following, we first derive the photon self-energies 
and then solve the resulting Dyson equation in the steady 
state. This corresponds to treating the electron-photon 
coupling in the random-phase approximation (RPA). It 
provides us with explicit expressions for the photon prop¬ 
agators - both for normal conducting and for supercon¬ 
ducting leads. Details of the calculation are shifted to 
Appendices A, B and C. 


A. Electronic feedback 




FIG. 2. Feynman graphs of the one-loop bosonic self-energy 
due to the coupling to conduction (c) and valence (v) 
electrons, see Eq. (30). External photon propagators (wiggly 
lines) are not part of the self-energy. Each vertex is associ¬ 
ated with a coupling constant |go|. Wiggly lines denote pho¬ 
tons, solid (dashed) lines denote conduction (valence) electron 
propagators. The anomalous contributions 1111,112*2 appear 
only for superconducting leads A^, A„ ^ 0. 


the transition of electrons between conduction and va¬ 
lence bands. 

The electronic part of the photon self-energy has the 
following structure in the particle-hole basis of $q: 


ri®q'(t,t') = 


e*‘^(*+*')nf\q(T) 

ni'i,q(rj 


nf2,q(r) 

e-*'^(‘+‘')n®2*2^g(r) 


(28) 


where r = t — t' denotes the time difference. This 
structure is identical for the retarded, advanced and the 
Keldysh components of the self-energy. The anomalous 
components on the diagonal, which are associated with 
bqb-q and b]jb^_^, depend on the absolute time t + t' via 
the phase 

(f>{t 1') = eVo{t 1') — -|- (/>„ -|- 2(j)g . (29) 

Here, (j)c/v = arg(Ac/„) denote the (constant) phases of 
the superconducting gaps = |Ao,|e®‘^“, which depend 
on microscopic details at the initial time t = 0, and (pg 
is the phase of the coupling constant go = \go\D^° ■ The 
remaining parts nf\ ^ and n |*2 q, as well as the normal 
components on the off-diagonals nf* 2 _q,n|\ q, only de¬ 
pend on the time difference r = t — t'. They are also 
independent of the phases 4>c,v,g- 

We treat these self-energies in the one-loop approxima¬ 
tion. As shown in detail in Appendix A 1 the self-energies 
are given by 


Let us first focus on the electronic contribution to the 
effective action in Eq. (27). It is described by the self¬ 
energy ri®*(t,t') and contains the information about the 
photon absorption and emission processes that involve 


\^^XLi^:t')= (30) 

-z^tr[7“ 9tit)Go,kit,t')^ti g^{t')Go,k+qit',t)] . 

k 
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The corresponding Feynman diagrams are shown in 
Fig. 2. The remaining task is to insert the expression 
for the fermionic Green’s functions, perform the traces 
over the Keldysh and Nambu indices and to carry out 
the summation over electronic momenta k. 


1. Energy scales in the system 

Before we describe the explicit result for the self¬ 
energies in the next sections, we want to discuss the dif¬ 
ferent energy scales that are present in the problem. To 
give a numerical estimate of the different scales, we use 
realistic values of GaAs. 

There are five important energy scales in the problem: 
(i) the applied bias voltage eVb, which is of the order of 
the semiconducting bandgap D. This is the largest en¬ 
ergy scale in the problem. It sets the photon energy 
Hu!q « eVb and the photon momentum q = eVo/hc, 
where c denotes the speed of light. For GaAs this en¬ 
ergy scale is given by hujq = 1.424 eV; (ii) the semi¬ 
conductor Fermi energy ep. It is determined by the 
effective masses mc,TO„ and the carrier density 
as eF,a = (37r^)^/^Ii^na^^/(2TOQ). In GaAs one finds 
rric — O.OGTrUe and two hole bands with = 0.45me 
and my ^2 = 0.082me, where rrie is the bare electron 
mass. Typical values for the carrier density are n = 
10^®cm“^ — 10^®cm“^. For this range of doping and 
approximately using TOc = —= O.ObTwe, we find 
ep = 2.5 • 10“^ — 2.5 • 10“^ eV. The Fermi velocity follows 
to vp = \j2eplm = 4 • 10“"^c — 4 • 10“^c; (hi) the (prox¬ 
imity induced) superconducting gap |A| for which we es¬ 
timate |A| Ri ImeV; (iv) the strength of the electron- 
photon coupling IsoPp with conduction band density of 
states p. For the range of carrier densities above one finds 
the estimate = 8.3 • 10“®eV — 8.3 • 10“® eV^®. This 
is the smallest energy scale in the system, which justifies 
our perturbative approach; (v) the coupling to an exter¬ 
nal photon bath induces a photon decay rate t\ (for each 
photon mode q). Our assumption of a steady-state with 
a finite photon number requires that rj must be larger 
than the LED photon production rate. 

A result of this estimate is that one cannot neglect the 
photon momentum q. Although it is much smaller than 
typical electronic momenta |q|/|fcF| = eVocKvpmcC^) « 
0.01 — 0.1, its associated energy scale is of the same order 
as the superconducting gap 

UF|q|«|A|. (31) 

In the following, we discuss the explicit result for the elec¬ 
tronic part of the photon self-energy We first 

describe the case of normal conducting leads and then 
the case of superconducting leads. Taking a finite pho¬ 
ton momentum into account serves as a physical cutoff 
for divergences that would otherwise occur in the imagi¬ 
nary part for the superconducting case. 



FIG. 3. (color online). Schematic of the symmetric electronic 
band dispersion model for normal conducting (left) and su¬ 
perconducting leads (right). Normal state conduction and 
valence band dispersion obey ^c(fe) = —^v{k), vp denotes the 
Fermi velocity, kp the Fermi momentum, B the filling factor 
of the bands and eFo is the applied bias voltage. Electronic 
transitions involving photon emission (absorption) are possi¬ 
ble at photon energies ojq < eVo {ujq > eFo)- Photon momen¬ 
tum q is properly taken into account and results in electronic 
transitions that are not vertical. In the presence of supercon¬ 
ductivity the electrons at the Fermi energies are gapped out, 
resulting in allowed transitions for \ujq — eVb| > 2A. 


(arbitrary units) 



FIG. 4. (color online). Normal retarded photon self-energy 
112 * 1 ’^ (a;) at fixed momentum q as a function of frequency 
LJ. For normal conducting leads the imaginary part is shown 
in dashed green. The self-energy changes sign at uj = eVb; 
positive (negative) sign corresponds to photon production 
(absorption). For superconducting leads the imaginary part 
(blue solid) exhibits jumps at a; = eVb ± 2|A| (here |Ac| = 
|A„| = |A|), remains constant over an interval of width 
{vp\q\‘^/2)^/\A\ and then falls off towards the normal state 
solution. For superconducting leads the real part (red dashed) 
exhibits logarithmic divergences at the position of the jumps 
in the imaginary part. 


8. Normal conducting leads 

Let us first discuss the electronic contribution to the 
photon self-energy in the case of normal con¬ 

ducting leads. In this case, the diagonal elements in 
Eq. (28), which are proportional to the product Ac Ay, 
vanish: = n| 2 ^q(t,t') = 0. 

To be able to perform the summation over momenta 
in the off-diagonal elements analytically, we make the 
symmetric choice of assuming parabolic conduction and 
valence bands bands with effective masses uic = — 
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and equal superconducting gap amplitudes |Ac| = |A„| = 
|A|. As a result, the electronic dispersions fulfill = 
—^y{k). The fermionic density of states (DOS) Pc(w) = 
Pv{uj) = p{uj) is only weakly energy dependent around 
the Fermi energy, where p(0) = pp corresponds to the 
DOS at Pc,v This situation is schematically depicted 
in Fig. 3. For general dispersions the summation can 
easily be performed numerically. As shown in detail in 
Appendix B, one finds for the retarded and advanced 
self-energies at Tp = 0 


n 


cl,fl 

21.q 


(w) 



'VF\q\ 

sign(w_) 


for jw-l < VF\q\ 
for |a;_| > VF\q\ 




(32) 


Here, the frequencies a;± = wieVb are measured relative 
to the applied voltage. The retarded function is shown in 
Fig. 4. Finite fermionic temperatures Tp > 0 will smear 
the zero temperature results, but do not yield qualita¬ 
tively different results. As the fermions are assumed to 
be in equilibrium at temperature Tp, the corresponding 
Keldysh self-energies are given by 


nf2;^(w) = coth(^) 


nf2^(a.)-nf4(a;) 


(33) 


As can be seen easily the Keldysh self-energies obey the 
symmetry n° 2 ^(a;) = Eq. (32) we have 

neglected the real part which depends only weakly on 
frequency. It gives an unimportant renormalization of 
the photon resonance frequency ojq. The imaginary part 
of the self-energy ^(a;) describes the production and 
decay of photons due to the coupling to the electrons. 
From Fig. 3 we see that for T = 0 there are two possible 
transitions: (i) Photons with energy Uq < eVo are emit¬ 
ted due to transitions of an electron in the conduction 
band to the valence band, or (ii) photons with energy 
ujq > eFo are absorbed by raising an electron from the 
valence band to the conduction band. The absorption is 
associated with a negative imaginary part of the retarded 
self-energy in Eq. (32), while the emission is associated 
with a positive imaginary part. The linear dependence 
of n 2 *h^(w) around w = eVg arises due to the restricted 
phase space of decay and absorption processes for pho¬ 
tons with energy \ujq — eVbl < VF\q\- For larger energies 
\ujq — eVo\ > VF\q\ on the other hand, the self-energies 
reach a constant value because we linearize around the 
Fermi energies. 

We have seen that the production of photons is de¬ 
scribed by a retarded self-energy n|\ ^(w) with a pos¬ 
itive imaginary part for frequencies 0 < a; < eVg. It 
changes sign at the applied bias voltage eVg- If there was 
no further contribution to the total photon self-energy 
this would be nonphysical, since it does not correspond 
to the analytic structure required by a retarded bosonic 


self-energy, which must have a negative imaginary part 



FIG. 5. (color online). Anomalous retarded photon self¬ 
energy nf]’^(tj) as a function of frequency oj. Real part (red 
dashed) and imaginary part (blue solid) are non-zero only for 
superconducting leads and exhibit features similar to the nor¬ 
mal self-energy (see Fig. 4). Note that features occur around 
u) = 0 here. 


for w > 0. Physically, this corresponds to the fact that 
the assumption of a steady-state breaks down if photons 
are produced at a constant rate but there is no photon 
decay mechanism considered. The photon number would 
simply grow without bound. Therefore, in Sec. Ill B we 
take a coupling of the photons to an external photon 
bath into account. The associated bath induced photon 
self-energy n)’j)q*'(a;) will be added to the electronic con¬ 
tribution to the self-energy. The total self-energy will 
then have a negative imaginary part at positive frequen¬ 
cies w > 0 as required. 


3. Superconducting leads 


In the superconducting state the anomalous diagonal 
elements of the photon self-energy in Eq. (28) are non¬ 
zero. They depend on the total time T = (t + t')/2 and 
it is thus convenient to perform a Wigner transformation 
/(w, T) = fZ, dTf{T, T)e*-^ which yields 


n 


ol.fl 

q 


{uj,T) 


(34) 


with phase </>(2T) = 2eVgT + (j)y — (j)c + 2(j)g. The self¬ 
energy elements obey the relation n)’\’'^(a;) = 1422^(‘^) 
and nf4(a;) = n“\’;^^(-a;). 

As shown in Appendix B 2, the one-loop self-energies 
at Tp = 0 take the form 












(35) 


^cl,R/A 

^^11,q 


j-rCl,R/A 

^hl,q 



Rk,v'^k,vRkAq,c'^kAq^c 

Rk,v'^k,vRk-\-q,c^k-^q,c 

1^012^ 
k 

Lu — Ey{k) — Ec{k + q) + i0 

UJ + Ey{k) + Ec(k + q) ±i0 






1^012^ 
k 

UJ- 

— Ey{k) — Ec{k -|- q) ± iO 

UJ- + Ey{k) + Ec{k + q) ±i0 


(36) 


where Ufc,„ = and Vk,a = - 4^) 

are the superconducting coherence factors for the con¬ 
duction and valence band. They contain the Bo- 
goliubov quasi-particle dispersion relation Ea{k) = 
v/?„(fc)2 + |A„|2. 

The normal component is shown in Fig. 4 and the 
anomalous component in Fig. 5. From the retarded self¬ 
energies Hi]’g (w) we easily get the Keldysh self-energies 
as 


n 


cl,if 
11.q 


(c.) = n 


c\,K 

22,q 


cl, 

-21,g 


(ai) = n 


c\,K 

12,q 


(uj) = coth(^) [nii^(a;) - nii;^(a;)] , 

(-w) = coth(^) [n2i^(w) - n2‘i;^(a;)] 

(37) 


Let us discuss the retarded self-energies in some de¬ 
tail. The combination Uk^a'Vk,a = \^a\/{‘2‘Ea{k)) and 
the diagonal elements sjce therefore proportional 

to the product |A„||Ac|, i.e., they are non-zero only 
if both leads exhibit superconductivity. In the follow¬ 
ing, we assume for convenience a momentum indepen¬ 
dent gap function which is identical for the two bands 
I Ac I = |A„| = |A|. Note that the phases of the super¬ 
conducting order parameters ipc and 4>y are factored out 
explicitly in Eq. (34). 

Both diagonal and off-diagonal elements of the self¬ 
energy show similar behavior. Their imaginary part van¬ 
ishes in a region of width 4|A|. While the normal compo¬ 
nents are zero for |a;_| < 2|A| the anomalous components 
vanish for |a;| < 2|A|. At the border of these regions the 
functions exhibit a jump in the imaginary part, then they 
remain constant over a frequency window of the order of 
i;|,|qp/|A| and finally decay towards the result for nor¬ 
mal conducting leads (see Eq. (32)) further away from 
the gapped region. 

Explicitly, the real (imaginary) part of the normal self¬ 
energy 1121 (anti)symmetric around the region 

u! = eVg. The imaginary part is zero for jw —eVbl < 2|A|. 
At the border it exhibits a jump of size 


lim 
(5—)-0 + 


Imn^\’«(eyo + 2|A|+(5) 


2 


Iffol^ Pf 


|A| 

VF\q\ ' 
(38) 


Here and in the following we consider ?;i?|q|/2 < |A|, 
which holds for our realistic choice of parameters (see 
Sec. IIIAl). A jump in the imaginary part yields via 


the Kramers-Kronig relations a logarithmic divergence 
in the real part at w = eVo- If we define the functions 


R{w,x) = 


2 X 


11 + 


X 


\9of —^i\w\ - 2) arcsin 
2 X 


2-|u;|| y|2- \w 

min[x, ^/\w\^ 


V\M - 2 


(39) 


the self-energy for frequencies w « eVg ± 2|A| can be 
expressed as 


-cl. 

21.g 



w- VF\q\\ 
lAI’ 2|A| ) 


■ i sign(a;_)/ ^ 


w- VF\q\\ 
lAI’ 2|A| ) 
(40) 


Since \ggf Pf ^ |A| (see Sec. IIIAl), the real part is 
only important around its logarithmic divergence. It 
quickly decays towards the normal conducting result 
away from the resonance. The imaginary part of H^^ is 
given by Eq. (40) close to the jump and also approaches 
the constant result that we found for normal conducting 
leads (see Eq. (32)). 

The vanishing of the imaginary part for |w_| < 2|A| 
can easily be understood if we look at Fig. 1. For zero 
temperature there is no possibility for one photon to 
be either absorbed or emitted within this energy range 
{eVo — 2|A|,eVb + 2|A|} due to the superconducting gaps 
in the conduction and the valence bands. On the other 
hand, in a superconductor the DOS diverges at energies 
±|A| relative to the Fermi energy and there are thus 
many electronic states leading to an enhanced emission 
and absorption of photons coupling to those states. The 
imaginary part of the self-energy is therefore enhanced in 
this region compared to the normal conductor. 

We note that it is absolutely essential that we take the 
finite photon momentum q into account. Otherwise, the 
imaginary part would exhibit a square-root divergence 
at |a;_| = 2|A|. This divergence is cut-off by finite q at 

||a;_| — 2|A|| = leading to a finite jump instead. 

Since the imaginary part corresponds to the photon pro¬ 
duction rate in the superconducting LED, this rate would 
diverge if one neglects the photon momentum. 

The anomalous diagonal components of the self-energy 
nf- are shown in Fig. 5. Close to |a;| = 2|A| they can 
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also be expressed by the functions defined in Eqs.(39) as 


with component 


el, 

11,q 



UJ 

|A|’ 2|A| ) 


+ i sign 



oj 

|A|’ 2|A| ^ 


(41) 


Away from |a;| « 2|A| the imaginary part decays like 
1/w^ to zero. This is a faster decay than predicted by 
Eq. (41). 

At finite but small temperatures Tp ^ |A|, these re¬ 
sults acquire small corrections such as as an exponen¬ 
tially suppressed imaginary part in the gapped regions. 
At higher temperatures the described features are sup¬ 
pressed, but this also leads to the breakdown of super¬ 
conductivity and one approaches the results for normal 
conducting leads. 


B. Influence of the photon bath 


To achieve a steady-state in the system it is required 
that the photons, which are produced by electrons mak¬ 
ing a transition from conduction to valence band, may 
also be absorbed. We thus consider the coupling to an ex¬ 
ternal photon bath as a decay mechanism^^’^^. As shown 
in Appendix A 2, the coupling to the bath (see Eq. (6)) 
gives rise to an additional contribution to the photon 
self-energy, which corresponds to a finite photon lifetime 
T~^ — r). If we neglect the unimportant real part, which 
has no divergent features, the bath induced photon self¬ 
energy reads 


l=lbath.fl(^) 


/ Q _nbath.F(_^)\ 


(42) 


nbath.iC(^) = coth(^) _ nb-“'’^(w)] 

(46) 

would diverge as l/w for small w. Here, we have assumed 
that the external photon bath is in thermal equilibrium 
with photon temperature Tg. One can easily incorpo¬ 
rate a different external photon distribution by replacing 
coth(a;/2TB) with an arbitrary bath distribution func¬ 
tion B{u!) = 1-1- 2n'^’^‘’'(a;), where n*’®'**'(w) denotes the 
number of bath photons in a state of energy w. 


C. Dressed Photon propagators 

To calculate observable quantities such as the lumi¬ 
nescence or the statistical properties of the light emit¬ 
ted from the superconducting LED, we need to find the 
dressed photon propagator Dq. The photon self-energies 
n®' and nbath determine the dressed photon propaga¬ 
tor Dq via the Dyson equation, which reads in Keldysh 
space^® 

{Do,l-^^)oDq = i, (47) 

where Hq = Hg -|- 11®*^®®. Here, o denotes a convolution 

in time. Explicitly, this corresponds to three coupled in¬ 
tegral equations for the retarded, advanced and Keldysh 
components of the dressed propagator 

= i (48) 

([I)o^q]-i-n^)ol)f =fl^oD^. (49) 


with 

nbath.fl(^) ^ ^ ( 43 ) 

We assumed that the coupling App/ = A(wp) depends 
only on the photon frequency and introduced the DOS of 
the external photon bath /5bath(<j-'). The imaginary part is 
determined by the spectral function of the external bath, 
which we assume to be of Ohmic form 


One can omit the infinitesimal component ^ fO, 

which acts as a regularization in the non-interacting sys¬ 
tem, since the coupling to the electrons and the bath 
induces a finite Keldysh self-energy. Since the anoma¬ 
lous components of the self-energy n®^(T, T) depend on 
the absolute time T, we perform a Wigner transformation 
with the ansatz for the full propagator (see Appendix C 
for details) 


7r|A(a;)|Vbath(w) = ry0(w)^^-^ (44) 

with A ^ eVb- The spectral function is constant rj > 0 
for w ^ A and a vanishing function for a; —>■ 0. The fact 
that it decays to zero at small frequency is important 
since otherwise the Keldysh self-energy 


H 


bath./C 


(w,r) 


/ 0 n®^t®’^(-a;)\ 

(^nbath,F(^) 0 ) 


(45) 




(50) 


with the phase factor given in Eq. (29). The dependence 
on the absolute time is similar but of opposite sign to the 
one found in the self-energy matrices. The explicit forms 
of the full retarded and advanced photon propagators 
read 
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— ^22,qi^) — 




[w + eVo -UJq- ll 2 i^q 




-Wq-n« q(a; + eVb)][a;-eyo+Wq . 

(jj — 2eVo + u)q + — ^eVg) 


+ nf2,^(u; -eyo)] + [nfi,^M]' 


_T_ 

[uJ-UJq- nfi q(w)] [w - 2eVo +UJq+ nf 2 ,q(w - 26^0)] + \tif^l^q{uJ - V)] ' 


The Keldysh components are given by 

^n.qM = -[D^2,qi^W = D^^qi^) [flfi,qM,q(o;) + nf2,q(w - V)D^,Ju; - F)] 

+ i 5 f 2 .q(^ + [^21.qi^ + V)DA.qi^) + ^22 .qi^) K .qi^ ” ^)] 

Dl2.,q{^) = -[D21..q{-0o)r = D^l.qioJ - V) [nfi,q(cC - V)Dt2,q{uj) + U^^.qi^ - 2V)b^^^q{u: - F)] 

+ i5f2,q(w) Mi,<ii^)Dt2,q{^) + flgju; - F)44,q(a; - F)] 


(51) 

(52) 


(53) 

(54) 


By inverting the Dyson equation we have summed up 
the complete RPA series of bubble diagrams with the 
electronic contribution to the photon self-energy given 
in Fig. 2 and the bath contribution given by Eq. (43). 
In fact, the expressions for the dressed propagators are 
formally exact if the self-energy was known exactly. This 
follows from the fact that the structure of the self-energy 
in Eq. (34) holds to all orders in perturbation theory.. 
By inverting the Dyson equation we have only employed 
this general structure. 

From Eq. (52) we find that the photon particle prop¬ 
agator D ^2 q ('^) shows features at w « Wg and to « 
eFo ± 2|A|. In the case that the photon is off-resonant 
with the applied voltage, which is determined by the 
semiconductor bandgap, ojq ^ eF) ± 2| A| the photon can 
just propagate through the system and is only weakly 
interacting with the electrons. In contrast, if the pho¬ 
ton is resonant Wq « eFo ± 2|A|, the enhanced spectral 
weights around ui « eFo i 2|A| correspond to photon- 
exciton bound states or polaritons. Here, a photon can 
excite a quasi-particle from the superconducting edges of 
the valence band to the conduction band, which again re¬ 
combines under emission of a photon. This process may 
repeat itself an arbitrary number of times thus forming 
an electron-photon bound state, a polariton. Since the 
DOS of the two superconducting bands diverge at ±|A| 
(measured from the the Fermi energy), only scattering 
processes with a photon matching the energy difference 
Wq « eFo ± 2|A| give rise to a large effective coupling 
between photons and electrons and to the formation of 
polaritons. 


IV. LUMINESCENCE AND SQUEEZING 
PROPERTIES OF THE SUPERCONDUCTING 
LED 

In this section, we investigate the photon luminescence 
C[(jjq) = (6qfoq) of the superconductiug LED, z.e, the 
number of photons present in the system in the steady 


state. We consider both the case of normal conducting 
leads and the one of superconducting leads. Supercon¬ 
ductivity leads to a strong enhancement of the lumines¬ 
cence in a frequency window close to Wq = cFq — 2|A|. 
We demonstrate that the superconducting LED emits 
entangled photon pairs and produces squeezed light of 
frequency ujq = eFo. The squeezing occurs in certain 
two-mode quadrature operators of the light field defined 
below and implies that the fluctuations in one of the 
quadrature components falls below the minimal uncer¬ 
tainty of coherent states. 


A. Photon luminescence 

The photon luminescence is defined as the expectation 
value ( bqbq), which can be expressed via the lesser pho¬ 
ton Green’s function Db,qi^b') = ~'>-{bq {t)b~ {t')) as 

/ OO 7 

^T)<2_q(w) (55) 

The lesser propagator D^.q = U^i 2 ,q “ Di 2 ,q + Di 2 ,q) 
is a linear combination of the retarded, advanced and 
Keldysh propagators given in Eqs. (51)-(53). 

In Fig. 6 we present the luminescence in the normal 
and in the superconducting state for an external photon 
bath at Tg = 0. In Fig. 7 we show the luminescence in 
the presence of bath photons for a bath kept at temper¬ 
ature Tb = eFo/2. Those photons can be absorbed in 
the semiconductor junction and transfer electrons from 
valence to conduction bands. 

Let us first focus the case of Tg = 0. For normal con¬ 
ducting leads we then observe that only photons with 
frequency Wq < eVo are present in the system. This fol¬ 
lows from the fact that electronic transitions are only 
available with energy difference ujq < eFo (see Fig. 3). 
The frequency scale on which the luminescence increases 
from zero to its constant value is given by the photon 
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FIG. 6. (color online). Luminescence C{u}q) = (bqbg) as a 
function of photon frequency cUq = c|q|. We show the lu¬ 
minescence of the LED with normal conducting leads (blue) 
and superconducting leads (red). Temperature of the exter¬ 
nal photon bath is kept at Tb = 0. Enhanced luminescence 
around ojq = eVb — 2|A| in presence of superconductivity is 
due to quasi-particles tunneling from conduction to valence 
band and larger density of states at the edges of the super¬ 
conducting gap (see Fig. 1). Inset enlarges region around res¬ 
onance ujq = eVb showing luminescence peak due to Cooper 
pair tunneling. 



-1-^- I ' I ^ 

eVo - 2|A| eVo eVo + 2|A| <^q 


FIG. 7. (color online). Luminescence C{ujq) = {b\bq) as a 
function of photon frequency ujq = c|q|. External photon bath 
is kept at temperature Tb = eVb/2. Thermal bath photons 
(green dashed) with frequency ujq > eVb -I- 2|A| can be ab¬ 
sorbed by the p-n junction both for normal conducting (blue) 
and for superconducting leads (red). Inset enlarges frequency 
region where the Cooper pair peak and the absorption can be 
seen. 


momentum vf\q\ with |q| = eVo/c. In the case of super¬ 
conducting leads, we clearly observe a strong enhance¬ 
ment of the number of photons with frequency close to 
uiq = eVg — 2|A|. This derives from the large number 
of electronic states that are pushed out of the gapped 
region in the superconductor to the border. Photons 
are produced via a recombination of (Bogoliubov) quasi¬ 
particles from the conduction band to the valence band as 
depicted in Fig. 1. The behavior of the luminescence can 


be traced back to the photon self-energy shown in Fig. 4, 
if one notices that the imaginary parts of the self-energies 
correspond to photon production rate (for positive imag¬ 
inary part) and decay rate (for negative imaginary part). 

In addition to the enhanced photon production at the 
band edges, the luminescence also exhibits a Cooper pair 
peak at coq = eVg in the presence of superconductivity. 
This stems from Cooper pairs that are transferred from 
the conduction to the valence band. This process is of 
the order |qo|^ and arises from the diagram shown in 
Fig. 8 that contains two anomalous self-energy bubbles. 
This contributions is taken into account in the RPA sum¬ 
mation. The luminescence peaks are characterized by a 
width Su} = max[77, |A|] which is also the width 

of the plateau of the imaginary part of the retarded pho¬ 
ton self-energies in Eq. (36). 

To obtain the numerical result in Fig. 6, we have used 
parameters that are consistent with our general discus¬ 
sion of energy scales in Sec. Ill A1. Specifically, we have 
expressed all energies in units of the superconducting 
gap |A| (a realistic value is |A| = ImeV). We have set 
\go\^ PF = |A|/50, vf/c = 0.001, and ePg = 1000|A| = 
leV. We further assume a simple quadratic conduction 
band with a band edge that lies at a distance B = Vg/10 
below the chemical potential The electronic density 
of states thus reads p{e) = pf\/^ + For a non-lasing 
steady-state to exist, the photon decay rate due to the 
bath f] must be larger than (see Eq. (38)) 


= lgol%f |A| 
2{vF/c)ujq 


(56) 


We choose p = 1.5?7min as even larger decay rates (in 
particular p > |A|) result in a photon linewidth that is 
larger than the superconducting gap which smears the 
features observed in Fig. 6. As soon as p < p^in the 
system exhibits lasing and the photon number diverges 
at a particular frequency. We discuss this possibility in 
detail in Sec. V. 

If the bath contains real photons they can be absorbed 
by the hetero-structure. This case is shown in Fig. 7 
where we assume a thermal photon bath at temperature 
Tb = eVg/2. Photons of frequency ojq > eVg + 2|A| may 
be absorbed promoting electrons from the valence to the 
conduction band. This occurs both for normal conduct¬ 
ing as well as for superconducting leads. Here, we neglect 
the effect on the fermionic distribution and assume that 
the fermionic distribution functions remains at Tp = 0. 
This can be realized in practice by a fermionic bath of this 
temperature. In addition to the photon emission peaks, 
which are unaffected by the presence of bath photons, 
we now clearly observe an additional absorption dip at 
ujq = eVg + 2|A| in the presence of superconductivity. 
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FIG. 8. Feynman graphs for photon emission processes of 
order CI((?o) that contribute to the luminescence C{ujg). Part 
(b) contains two anomalous self-energy loops 1111,1122 and 
gives rise to the Cooper pair peak seen at ojq = eVb- 



B. Light Squeezing 

In this section, we investigate the statistical proper¬ 
ties of the light that is emitted from the superconduct¬ 
ing LED. We find that the photons inherit the anoma¬ 
lous correlations that are present between electrons in a 
superconductor. In particular, the photon expectation 
value {b^jb'Lg) is non-zero, which reduces the quantum 
fluctuations of certain two-mode quadrature operators of 
the light field to a value below that found for a coher¬ 
ent state^®. The fluctuations of the conjugate operator 
are of course increased such that the Heisenberg uncer¬ 
tainty limit is obeyed. The superconducting LED there¬ 
fore emits two-mode squeezed light, which can act as a 
resource for quantum information processing and metrol- 

The light squeezing in this setup was studied in Ref. 16 
within lowest order perturbation theory in the electron- 
photon coupling go. To this order, the luminescence 
{bjjbq) vanishes on resonance ojq = eVo (at Tb = 0) 
where the lowest order terms are O{go). The amount 
of squeezing is then solely determined by the anomalous 
expectation value Here, we consider the effect of 

the diagonal (photon number) expectation values {bjjbq) 
on the squeezing by summing the complete RPA series. 
We make a quantitative prediction for the maximal re¬ 
duction of quantum fluctuations that is achievable in this 
setup. 

Squeezing occurs in the two-mode quadrature opera¬ 
tors 


FIG. 9. (color online). Anomalous photon expectation value 
CAiijJq) = e’'’^^^^\bq{t)b-q{t)) US u function of photon fre¬ 
quency iJq = c\q\. We show the real part (red), imaginary 
part (blue) and absolute value (green). External bath tem¬ 
perature is set to Tb = 0. 

below the Heisenberg uncertainty limit: ((AA^)^ ) ^ 3 
for either i = + or i = —. If the number of photons 
( bjjbq ) in a mode q is zero, a finite anomalous expec¬ 
tation value of Re ( bqb-q ) always results in squeezing. 
On the other hand, in the presence of photons one must 
compare the expectation values {bjjbq) and (feg^^Lg) to 
each other. Since the number of photons is smaller at 
low temperature, the resulting squeezing amplitudes are 
larger. We thus focus on the case of Tb = 0 in the fol¬ 
lowing. A crucial observation is that close to resonance 
Wg = eVo it holds that the luminescence ( bjjbq) = O{go) 
while the anomalous expectation values are already of or¬ 
der ( bjjbj^^ ) = O{go). We obtain the anomalous photon 
expectation value 

(6g(t)Lg(f)) ^^D<JuJ,T = t) 

= e*[2“A-0(2t)] f (59) 

J 27r 

by an integral over the lesser anomalous propagator 
~ ~ ^ii.q + ^^.q)' Wo define the anoma¬ 
lous luminescence by CA{uJq) = (6g(t)6_g(t)) = 

/ ^iDii q{uj) = e®*‘'|£A(wg)| to arrive at 


A^=M±[bl + blq±h.c.], (57) 

where Af+ = 2“^/^ and Af- = The photon 

operators are defined in the rotating frame bq = bqe^‘^A^ 
The fluctuations of these two-mode quadrature operators 
are given by^® 


(AA 




1 -b 2( bjjbq ) ± 2Re ( bqb-q 


(58) 


where (AA)^ = (A — (A))^. The state is truly squeezed 
if the fluctuations in one of the quadrature operators fall 


((AA±)") = l(l + 2£(a;g) (60) 

± 2cos[2(wg - eVo)t + (j)q] |£A(wq)|) • 

Here, the initial phase of the last terms reads (j)q = 
arg[A„A*(/Q£^(a;g)] = (j)y — (j)c + 2(j)g -b Iq and the time 
dependence vanishes for photons on resonance Wg = eVb- 
The anomalous luminescence CA{<AJq) is shown in Fig. 9. 
At Tb = 0 it exhibits two main peaks: one at Wg = 
eVg ~ 2|A| corresponding to the transition of a Bogoli- 
ubov quasi-particle from the conduction to the valence 
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band and one at ujq = eVo corresponding to the transi¬ 
tion of a Cooper pair. It is important to note that the 
peak at ujq = cVq — 2|A| implies that breaking up (two) 
different Cooper pairs still leads to the emission of corre¬ 
lated and phase coherent photons due to the presence of 
macroscopic electronic BCS condensates. Both processes 
are depicted in Fig. 1. Due to the finite number of pho¬ 
tons at ujq = eVb + 2|A| which are produced in the system 
(see luminescence in Fig. 6) we observe a small peak of 
CAi^q) at Wq = eVb + 2|A|. This peak corresponds to an 
absorption process of a Bogoliubov quasi-particle from 
the valence to the conduction band. 

In order to obtain squeezing for the mode with photon 
momentum q, it is required that |£A(wq)| > |>C(a;q)|. 
Comparing the normal luminescence C{u)q) in Fig. 6 
with the anomalous luminescence CA(uJq) in Fig. 9, we 
find that squeezing is maximal for photons on resonance 
ujq = eVoi which corresponds to the transitions involving 
Cooper pairs. 

In Fig. 10 we show the amount of squeezing on res¬ 
onance at Tb = 0. It shows that squeezing can be 
controlled by the relative phase between the two super¬ 
conductors (j)c — The fluctuations of either Akl+ or 
Ayl“ fall below the Heisenberg uncertainty limit for a 
broad range of relative phases. The maximal amount of 
squeezing for our realistic choice of parameters is about 
10 percent. 

In the experimental setup of the p-n junction the rela¬ 
tive phase A(() = 4>c — 4>v depends on microscopic details 
such as the initial switch-on time. In an experiment it 
will be random from experiment to experiment but fixed 
within one run. If one places two superconducting p-n 
junctions in parallel in a SQUID geometry one can con¬ 
trol the relative phase 4>c — </>„ between the two junctions. 
Since this will essentially change the individual relative 
phases A())^, A(/)^ as well, this provides a way to change 
(pc — 4>v with a magnetic field. 


V. STEADY STATE AND LASING THRESHOLD 

In our calculation so far we have assumed that the sys¬ 
tem reaches a steady-state with a finite number of pho¬ 
tons in the system. This requires that the bath absorbs 
the photons that are produced in the LED sufficiently 
fast. In this section, we derive the exact requirements 
that have to be fulfilled for this to be the case. Other¬ 
wise, the luminescence diverges at certain frequencies and 
the system exhibits lasing. We discuss the lasing condi¬ 
tions and give the frequency window where lasing occurs. 
Finally, we calculate the steady-state photon distribution 
using rate equations. 


A. Lasing condition 

To derive the conditions that are required for a steady- 
state with a finite photon number to exist, we follow a 



FIG. 10. (color online). Fluctuations {(AAg)^) of the 
two-mode quadrature operators App for \q\ = eVofc (on res¬ 
onance) as a function of superconducting phase difference 
4>c — 4>v ‘itpg. The phase cpg denotes the constant phase 
of the electron-photon coupling constant go- For a range of 
phase differences one of the fluctuation amplitudes falls below 
the Heisenberg uncertainty minimum for the symmetric case 
{ (AAq ) ) = 1/4. This shows that a superconducting LED 
emits squeezed light. Inset shows mean and uncertainty of 
two-mode quadrature operators Ag . The orientation of the 
squeezing ellipse is controlled by pc — Pv 


discussion given in Ref. 41. The key idea is to deter¬ 
mine when and for which frequencies the luminescence 
exhibits a divergence. This corresponds to a transition 
into a lasing regime. We first consider the case of normal 
conducting leads and then the one of superconducting 
leads. 

The luminescence C{ujq) = (b^qbq) can be obtained via 
the lesser Green’s function 0^2,q = 5 (^i 2 .q ~ -^ 12 .q + 
^ 12 ,q) (see Eq. (55)). We can parametrize the inverse 
retarded, advanced and Keldysh Green’s functions as 


[^f2(q(^)]-'=^,(^)±*^.(^) (61) 

[l^f2.q(w)]-'=*C'q(a;). (62) 


Using the matrix structure in Keldysh space 
Eq. (20)) one obtains immediately 


= [Aq{u:)±iBq{u))] 

l^f2,q(w) = - 


[D^2,q{^)]- 


\D^2,q{^)]-^[Dt2,q{^)] 


-1 


(see 


(63) 

(64) 


The zeros of Aq{uj) describe the excitations of the sys¬ 
tem and Bq(u}) their linewidth. Let us assume that there 
exists a resonance at the renormalized photon frequency 
uj = oj* such that Aq{uj*) = 0. It is then required that the 
imaginary part fulfills Bq{uj*) > 0 to obtain a proper re¬ 
tarded Green’s function with the poles lying in the lower 
complex frequency plane. 

In Eqs. (51)-(54) we have given the propagators within 
the RPA approximation. For normal conducting leads, 
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eVo - 2|A| 




FIG. 11. (color online). Lasing condition for normal conduct¬ 
ing leads. Upper panel shows nearly momentum independent 
imaginary part Bq{u}) ~ Bq^(ui) (green dashed) of inverse re¬ 
tarded propagator [Z)|| q(tj)]“^. Lasing occurs for momenta 
larger than |qc| ~ eVb/c. Corresponding real part Aq(uj) is 
shown for different momenta q around q^. The photon bath 
decay rate is chosen q < tt |qo|^ Pf to fulfill the laser thresh¬ 
old relation for frequencies Uq < eVo- Lasing occurs for mode 
with momentum q = qc where both real part Ag^(uj) (red) 
and imaginary part Bq^{oj) have simultaneous zeros. Lower 
panel shows the luminescence £{uJg) = (6q6q) as a function 
of photon energy ujq = c\q\. At cjq < Uq^ the luminescence 
diverges, which denotes a violation of our assumption of a 
(non-lasing) steady-state with finite photon number. In the 
lasing regime ojq < ujq^ (light blue region) the system pro¬ 
duces photons at a faster rate than the absorption q due to 
the bath occurs. 


the retarded function reduces to q(w)] ^ = uj — ujq — 
nfi, q{u}) and thus 

Aq(uj) = U! - UJ* (65) 

Bq{uj) = -Imn 2 \’^(w) - Imn 2 i*q’'^(a;). (66) 

Here, uj* = uiq + ReH^^.q denotes the renormalized pho¬ 
ton frequency, where the real part of the self-energy is 
only weakly frequency dependent. The imaginary part 
Im 021 ^( 0 ;) is given in Eq. (32) and Imn 2 i*q’'^(a;) in 
Eq. (43). The bath contribution to the self-energy ful¬ 
fills Imn 2 i*q’^(a;) « —q < 0 for all frequencies uj. In 

contrast, the electronic contribution Imn 2 '(^(w) changes 
sign at Cli = eVb and is positive for frequencies below 
resonance uj < eVb (see Fig. 4). This follows from the 
electronic population inversion and describes photon pro¬ 
duction via transitions of conduction electrons to the va¬ 


FIG. 12. (color online). Lasing condition for superconduct¬ 
ing leads. Upper panel shows nearly momentum independent 
Bq{uj) (red dashed) and real part Aq(u!) of inverse retarded 
propagator [D^ 2 ,q{uj)\~^ for different momenta q with energy 
close to uiq ~ eVb —2|A|. The photon decay rate q is chosen to 
fulfill the laser threshold relation in Eq. (70) in the supercon¬ 
ducting state. The lasing regime is bounded by momenta qi 
and q 2 which correspond to the photon modes where Aq{uj) 
and Bq{uj) have simultaneous zeros (for some frequency lS). 
Note the strong renormalization of the photon frequency close 
to a; = eVb ~ 2|A| which increases the lasing frequency win¬ 
dow shown in the lower panel (in light blue). Lower panel 
shows divergence of luminescence for photon modes qi and 
q 2 and lasing regime (light blue) in between. 


lence band (see Fig. 3). Depending on the size of the 
bath decay rate q the total imaginary part Bq(uj) may 
now be negative for frequencies uj < eVg. Photon excita¬ 
tions with a resonance energy uj* for which this is the case 
violate the analytical requirements of a retarded propa¬ 
gator. This indicates the breakdown of our assumption 
of a steady-state with a finite photon number. At this 
point, the system exhibits lasing^^. 

This breakdown can also be observed as a divergence of 
the luminescence C{ujq) = (5^6q) = ^ dujD f2_q{^), 

where 


D 


< 

12 ,q 


(w) 


1 Cqjuj) - 2Bq{uj) 

2 Hq(w)2 -1- Bq{ujY ' 


(67) 


If the bath decay rate q < max |Imn 2 \’'^(a;)| « tt \gg'f' pp, 
the imaginary part vanishes Bq{p,q) = 0 for some 
frequency p.q = eVg — qvF\<i\/{'^\9o^ Pf)- Around 
Pq one can linearize Bq(uj) « I3{uj — Pq) with /? = 
TT |qo|^ Pf/(i’f|q|)- If the excitation energy w* « pq is 
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close to this critical energy, we can approximate the lu¬ 
minescence as 

1 du; Cg{u:)-2Bg{u;) 

^ 2j_^27rAg{ojy + B^iojy 

_ 1 /•“ du; CgifXg) - 2Bg{tXq) 

2 2 -k (w - w*)2 -h /?2(a; - /iq)2 

= . ( 68 ) 

4/3|w*-/iq| 

We have used that the dominant part of the integral 
comes from the region around w « /ig,w*. We note that 
the renormalization of the photon frequency is a small 
effect for normal conducting leads ujq « w*. 

As shown in Fig. 11, the photon number diverges like 
C{u!*) ~ (w*—/iq)“^ from above, where the luminescence 
is finite. This behavior is reminiscent of Bose-Einstein 
condensation in an ideal gas. The divergence occurs for 
photons with momentum Qc where the zeros of the real 
part Aq{uj) and of the imaginary part Bq(ui) occur at the 
same frequency w*^. Neglecting any frequency renormal¬ 
ization LOq = LOq, WO find 


\(lc\ 


eVo 

c + — 

' TrlsoTPi? 


(69) 


Due to its population inversion the electronic system in¬ 
duces an effective chemical potential = c|qc| for the 
photons. Within our approach, we cannot calculate the 
luminescence for frequencies Wq < Wq^ since our assump¬ 
tion of a steady-state with a finite photon number breaks 
down. The LED produces more photons in these modes 
as the external bath can absorb. For these frequencies 
one needs to perform perturbation theory around a las¬ 
ing state instead. 

To summarize, the system exhibits lasing if the imag¬ 
inary part of the retarded self-energy becomes positive 
for positive frequencies. The luminescence diverges for 
that photon mode qc for which both real and imaginary 
part of the retarded photon propagator exhibits zeros 
Aq^{oj) = Bq^{oj) = 0 for at least one uj. 

Turning to the case of superconducting leads, we ob¬ 
serve that the electronic contribution to the self-energy 
is strongly increased around uj = eVo — 2|A| (see Fig. 4). 
The laser threshold thus corresponds to a larger value of 
the bath absorption rate 



FIG. 13. Feynman graphs contributing to the rate equation 
of the photon number (see Eq. (72)). Solid (dashed) lines 
denote conduction (valence) electron propagators, solid dou¬ 
ble wiggly lines denote dressed photon propagators given in 
Eqs. (51)- (54) and dashed wiggly lines denote bath photon 
propagators. The first three diagrams represent the contribu¬ 
tion to the rate from emission and absorption of photons in¬ 
volving electronic transitions in the the superconducting LED. 
The fourth graph describes the photon exchange with the ex¬ 
ternal bath. 


the simultaneous zeros in the upper panel of the figure in¬ 
dicate, superconductivity leads to lasing in the frequency 
range Wq^ < Wq < Wq^ « eVb — 2|A|. The luminescence 
clearly diverges as l/(a;q — coq^ ^) at the border of this 
region. Within this region one needs to perform pertur¬ 
bation theory around a lasing state, which is beyond the 
scope of this work. 

To conclude, superconductivity leads to a sharp in¬ 
crease of the photon production rate of the LED around 
the frequency ujq « eVb — 2|A|. This induces a sharp 
increase of the electroluminescence C{uiq) and a smaller 
laser threshold 


B. Lasing conditions from analyzing rate equations 

The lasing condition in Eq. (70) can also be derived 
using rate equations. This approach is appealing due to 
its simplicity and physical transparency. It begins with 
calculating the change of the number of bosons hf, = 
J^qb^bq, which is given by 

r = = (71) 

where the time-dependent operators nh{t),H{t) are de¬ 
fined in the Heisenberg picture. The expectation value 
(O) of an operator O is defined with respect to the full 
Hamiltonian H in Eq. (1). The commutator can be eval¬ 
uated to 


luol^ Pf ) normal conducting leads 

\9o? Pf 2 y^\q\ 1 superconducting leads ’ 

(70) 

where 7r|A|/(2u_F|q|) > 1 for our realistic choice of 
parameters (see Sec. HIAl). In Fig. 12 we show 
Aq{ui),Bq{ijj) and the luminescence C{uj*) for supercon¬ 
ducting leads for tt |go|^ Pf < V < ^ 
there is no lasing in the absence of superconductivity. As 


[nb{t),H{t)] = [git)vi^At)bl,_^{t)ck>.a{t) (72) 

k,k' .(T 

-f ^ ', \^p,qbl,{t)cLq{t) — Xq^pap{t)bq{t)\ . 

p,q 

The expectation value of this commutator can be conve¬ 
niently evaluated within perturbation theory in go using 
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the Keldysh formalism^^. The total photon production 
rate T = rsLED+hbath is composed of two contributions: 
first, Tsled which describes the emission and absorption 
of photons due to electronic transitions in the LED and 


second, Tbath, which characterizes the photon exchange 
with the external bath. An explicit calculation, which 
takes the diagrams shown in Fig. 13 into account, yields 


Tbath = -25] [ —(73) 
Tsled = 2^ / —Imn 2 i;^(w)Imi:)^^ 2 ,q(‘*^) ([^ + nq{uj)]9{eVQ - w) + nq{uj)6{uj - eVo) 

q J V 

_ 1 ^ /■ . (74) 

TT ^^ J TT ^ 

Here, nq{oj) denotes the photon distribution function in the (superconducting) LED system and n*’‘^‘''(a;) denotes the 
photon distribution in the external bath. For a thermal bath, is given by the Bose-Einstein distribution 

function. The coupling to the bath drives the photon distribution in the system Uq towards the external bath 
distribution nq{uS) —> In contrast, the LED produces photons via spontaneous as well as stimulated emission 

for Lu < eVo and absorbs them for uj > eVg. 

We evaluate the integrals in Eq. (73) and (74) to leading order in the electron-photon coupling go- To this order, 
we only get an on-shell contribution ImZl (^2 q(^) = — ujq), because higher order terms in Imli(\ q{u>) ~ 5 q are 

neglected as the self-energies are already of order ©(^g). As a result, the frequency integration enforces nq{u)) = n{wq) 
and the total photon production rate simplifies to 




Iinn2\’^(a;q) 


[l -f n{u}q)]6{eVo - lOq) + n(a;q)6>(wq - eEg' 


(75) 


In a steady state the photon emission and absorption are balanced and the total rate T = Tq is zero for all photon 

modes q. Employing that Imnq®'**'’'^(wq) < 0 and sign[lmn 2 \’^(a;q)] ^ sign(eVg — ujq) we can infer the photon 
distribution in the steady state as 


l{ujq) = 


ImT[2i’^(a;q) 


9{eVo - ujq) + |lmn^'^‘^’-«(wq)| n'=’^*^(wq 


-Imn: 


bath,i? 


(wq) -Imn2i’^(a;q) 


(76) 


For ujq > eVg the distribution function is positive and 
finite. For ujq < eVg, however, the imaginary part of 
the total self-energy Imn|\ q(a;q) = Imnq^‘*'’^(a;q) -I- 

ImH^f q(wq) in the denominator of Eq. (76) may become 
zero. At this point, the production rate of the LED is ex¬ 
actly equal in magnitude to the decay rate into the bath. 
This results in the divergence of the photon distribution 
for energies ujq for which Imn^^ q(a;q) = 0. It follows 
that the laser threshold in the (superconducting) LED is 
given by our previous result in Eq. (70). 


VI. CONCLUSIONS 

We have investigated the electroluminescence and pho¬ 
tonic properties of a forward biased p-n junction in prox¬ 
imity to superconducting leads. We have shown that su¬ 


perconductivity leads to a significant enhancement of the 
luminescence of the light-emitting diode in a frequency 
window of the order of (uFeVg)^/(c^|A|) « 10“^eU. This 
effect stems from the increased density of states at the 
edges of the superconducting gap in the electronic band- 
structure of the valence and conduction bands. The in¬ 
creased photon production rate also reduces the lasing 
threshold in the system. By summing the complete infi¬ 
nite order RPA perturbation series in the photon-electron 
coupling, we were able to show that an additional lumi¬ 
nescence peak occurs on resonance due to the tunnel¬ 
ing of Cooper pairs from the conduction to the valence 
band. In addition, such a superconducting light-emitting 
diode emits two-mode squeezed light. The squeezing an¬ 
gle is controlled by the superconducting phase difference 
(pc — 4>v between conduction and valence band. This 
proves that one may transfer the macroscopic coherence 
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of an electronic superconducting condensate to photon 
pairs and manipulate the photonic coherence electroni¬ 
cally. 


ACKNOWLEDGMENTS 

We acknowledge useful discussions with P. Baireuther, 
U. Briskot, B. Narozhny, M. Schiitt and J. Schmalian. 
The Young Investigator Group of P.P.O. received finan¬ 
cial support from the Concept for the Future of the Karl¬ 
sruhe Institute of Technology (KIT), Germany, within 
the framework of the German Excellence Initiative. 


Appendix A: Effective photon action 

The unperturbed fermionic and bosonic part of the 
action in Eq. (7) can be written in the basis of the R,A,K 
fields as 

/ OO 

(Al) 

1 

^ q 

where we defined the fermionic and bosonic spinors and 
matrix Green’s functions in Section II. The electron- 
photon interaction and photon-bath interaction can also 
be written in terms of these spinors. 


1. Integrating out the electrons 


Let us consider first the electron-photon coupling 


Sin 


[dt(g{t)bk-k'{t)ck,cr{t)vk\cr{t)+h.cA = V2 9i{t)^k-k'At)] 

fc.fc' dc \ / Jc 2 


(A3) 


where c&fe = ($fc,i,$fc, 2 )^ 


{bki b-k)'^ and we define the matrices in the extended Nambu-space as 


A 

9 lit) 


1 

7 ^ 


/ 0 0 0 0 \ 

git) 0 0 0 1 

0 0 0 -git) 

V 0 0 0 0 / 


A 


92it) 


1 

71 


/o git) 
0 0 
0 0 
Vo 0 


0 o\ 
0 01 
0 0 
-git) 0/ 


The Keldysh contour can now be expressed by two integrations on the real axis for the -I- and — fields 

Sint = ^ j dt^a - [«'fc(i)]7X] 72 5i(f)$fc_fc,^,(f)] ^l,it) 

/ OO ^ A 

0 '“ 9iit)<^k-~k'A*)] '^k'it) 

K,K- °=q.<=i 


Ku'it) 


(A4) 


(A5) 


and finally we have rewritten the action in the R,A,K basis using 7 '’* = 1 and 7 '^ = dx in Keldysh space. Here, we 
use the convention that we always hrst evaluate the Keldysh matrix structure and after that the structure of the 
fermionic Nambu-space A or bosonic particle hole space o. Now, we can integrate out the fermions easily as we have 
a quadratic action 


gd‘S'c-t-5,;+S'int] _ 


D[^, 4']e 


^ giSpfe-l-trln[l+GoF] ^ 

Expanding the trace log we hnd in leading order in g that 

°° / I '\n-l-l 1 

trln[l + GoV] = J] -tr[(GoR)"] = --tr[(GoH)'] + 


(A 6 ) 
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k.k' 

-I f dtdt'j2tr\GoMt',t)[T. 


k.q 


i==l,’2 


/3=q,cl 
J = l,2 


^ f dtdt'J2 J2 ®q.iW 

q i,j^i,2 ^ k ^ 

Q;,^ = q,cl ^ 






where we still have to evaluate the trace over the Keldysh and Nambu structure. We get 


1 /“^ 

5'pL,eff = 2 H y 






with 




n 






fif^{t,t') = 




J 21, q'' 

I d,cl 

11,0 

I d,cl 

21,^'-'^’ " J22,q^ 

9r{t)Go,k{t,t')^>^ 9j{t')Go^k+q[t'G)] 




(A7) 


(A8) 


(A9) 


2. Integrating out the photon bath 


Following the same lines as the previous calculation, we want to integrate out the bath photons in (6) to find the 
corresponding self-energy for the SLED photons, see also Ref. [31 and 32]. It is straightforward to show that the 
corresponding self-energy is then given by 




n 


bath,A// -f/AX 

frbath/. ./N _ / G,t)\ 

RA 1 - t')) 




lAjRIK 

^bath.p 


[dGY 


dlLt,q(A) = 


U! — LOq zt iO 


dhath,q{^) — ['7bath,q(b-') “ '^iXth.q (b-')] 


(AlO) 


The assumed (SLED photon) momentum-independent coupling Xpp> = X{u!p') then leads to the simple relation (we 
omit the unimportant since featureless real part) 


fl^ath,iiM(^) 

l=[bath./f(^) 


/ 0 ±ir]{—oj)'\ 

\^Ti??(w) 0 ) 

Ro(a;)[nb-‘hA(^)_nbath.A(^)] 


(All) 

(A12) 
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with 


77(0;) = Im 


^ |A(wp)|^d{^^th,p(w) =7r^|A(j/p)|^(5(a;-a;p) = 7r J deiVbath(e) |A(e)|^ (5(w - e) 


= TT |A(a;)|^iVbath(w) 


(A13) 


where A^bath = ~ ^p) DOS of the photon bath. We choose an Markovian (or Ohmic) bath, where the 

DOS and the coupling are constant in frequency for all w > A (where A ^ eVo) and thus give rise to a constant 
absorption rate of the bath for the photons produced by the SLED (with characteristic energy uiq « eVo)- We use the 
parametrization 


? 7 (w) = r] ■ 


■ A 2 


e{uj) 


(AM) 


with 77 ( 0 ) = 0 (which is important as explained in the main text). The derivation from the constant behavior for 
small w <C eVb does not change the results and features of the SLED presented in this paper as we are interested at 
photons with ujq ~ eVb- 

The total effective action then reads 






2 ^ 

q ' 


dtdt' (t) 


£)o-i(t,t')-n«i(t,t')-ff 


bath 

q 






(A15) 


Appendix B: Calculation of the self-energy fl®’ 

We already derived the expressions for the bosonic self-energy caused by the superconducting electrons in the system 
in (A9). Let us now first derive the causal structure of this self-energy by performing the trace over the Keldysh 
indices 


[n®']f 5 ,q(t,i') = = -f^tr[7‘i gj{t')Gk+q{t'.t)] 


(19) 


= -i trAtr^ 


‘“' I 0 ID' V 0 


g.(^t') ( Go^k+qi^' 


G(bfe+q(^^ 0 / 


= G^fc(t,t') G^fc+,(t',t)+ 9,{t) G^^k{t,t') 9jit') G^fc+q(t',t)] 

k 

[n®^])j,q(iD') = -^'^trA[9^it) Go,kit,t') 9jit') Go,k+qit',t)+ 9iit) Go,kit,t') 9jit') Go,k+qit',t)] 


A. 


[n®'] 4 q(^D') = -i^trA 


■ ■ ^.K u ,r ^ 


9^it) Glkit,t') 9jit')Glk+qit',t) 

A rA „ , A . , T A 


A , 


- 9rit)[Glkit,t')- Gikit,t')] 9jit')[Glk+qit',t)- Gik+qit',t)] 


(Bl) 

(B2) 


(B3) 


which are just the usual forms of the retarded, advanced and Keldysh self-energies as known from standard text books. 
The next step will be the evaluation of the Nambu trace. Let us first write define the absolute anomalous propagator 
(a = R,A,K) 


Po,k,v/cit, t') 
Po,k,v/cit, t ) 


.iR/A 

\0,k,v/c 


(w) 


= e-'^^'^\P\o,k,,,/Jt-t') 

= e'^^^^\P\o.k,./cit-t') 

I ^c/v I 

(W ± * 0)2 - ^„/c(fc )2 - I 


(B4) 


2 
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where we separated the phase of the superconducting gaps |A^/c|- Defining also the phase of the 

electron-photon coupling go = |go|; the absolute time T = {t + t')/2 and relative time t = t — t' we can evaluate 


rA A 


= |5o|^ [|P|o.fe,„ (r) (-r) + (r) (-r)] 

trA[5r(t) Go%(M') 5.(t') [Gg;:(M')GSf,.c(i',i) + GS;:(M')GSf,,.(t',t)] 

= Igol^ + Gi%{r)Gl,^lf^J-r)] 

trA[52(t) G?,fc(t,t') 5i(t') Gl^+^it\t)] = 

- 150^ [GiTv(^)Gi^l'Lci-^) + gS’.cWgS’+.,.(-^)] 

+ + Po,kA^P')Po,k + q,vi*' P)] 


= e 


A, 


trA[52(t) G^AtA) 92{t') GA+A^'G)] = - 


9it)9{t') 


= |^^|2 (-r) + \P\A,c (r) (-r)] 


with the rotating phase 


^{2T) = 2eVoT + [A.A:^"] = 2eVoT - + </>. + 2cj)g 


(B5) 


(B 6 ) 


From (A9) and (B5) we see that we indeed find the structure presented in Eq. (28) in the photon particle-hole space 
o. Performing the Wigner transformation f{oj,T) = W’{/(t, T)}^^ = dT/(r, r)e*“'^, which is just a Fourier 
transformation in relative time r, we have the identity 


/ OO 7 

^^A(a;i)P(a;i-cu-5) 


In total, we can write down the photon self-energy part from the coupling to the electrons as 

^gi0(2T)nol,flM^^^ YfhRJA^^-^ \ 




with 


fjcl.R/A, , . I 2 

nii.q (w) = ^ |go 


duji 

^ I ~2A 


n 


cLR/A 

12.q 

c\M/A 


' —OO 
/•OO 


doji 


' —OO 
pOO 


T-fCLH/A/ X -I 2|Y^ / 

n 2 i,q M = -*|ffol 2 ^/ ^ 

^ J-oo 

_ /•OO 7 

Pt2^q^A) = i\9l\Yl I 


duJi 


l^^loX (^i) \P\o,k+q,c Ai -A + \P\o.k,. Ai) \P\ti+,,c Ai - A 

“ ‘^+) + APP^'Ak+HlAl - ‘*^+) 
AAP^o',k+q,cA^ ~ + Gg^fc.f Al)GA]^Aq^Al - ^-) 

\P\elttv l^l^fc+q.c (‘‘^l - ‘^) + \P\e,k,v (^l) \P\o!t^q,c (‘^l “ A 


where we defined uj± = w ± eVo and used several symmetry properties 

^{p),R/A/K, , _ _^ih),A/RIK, , 

^0M,v/c '.‘^1 “ ^0,k,v/c 1 

iKZI (A = wCT. (-A 


i0,k,v/c 


(B7) 

(B 8 ) 

(B9) 

(BIO) 

(Bll) 

(B12) 

(B13) 

(B14) 


of the fermionic propagators. The corresponding Feynman graphs are presented in Fig. 2. Using the symmetries 
(B12) it is easy to show that 

nn’^^-^(a;)=n 22 ^/^(a;) (B15) 

(B16) 

The Keldysh self-energies can be easily shown to be given by the thermal equilibrium expressions (37). 
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FIG. 14. Linearization of dispersion at the Fermi surface and finite momentum transfer by photon. 


1. Normal conductor 


In the normal conductor the propagators simplify to \P\^fc — 0 and 
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(BIT) 


where npi^j) is the Fermi function of the system, such that we find the standard results for the occurring particle-hole 
bubbles (the anomalous self-energies Iln q = 1122 ,q vanish in the normal conductor) 
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As depicted in Fig. 3, let us now consider an isotropic electron band ^c(fc) = ■fcd^l) and a symmetric hole band 
^„(fc) = —^c(fc) and linearize around the Fermi surface 

^c(fc) « 0) ■ [k- fc_F((/), 6 »)] = VF{k - kF) (B19) 

where (/>, 9 parametrize the orbital direction of k in spherical coordinates and we can define the DOS as usual p(^) = 

<5[^ —Cc(fe)] ■ Let us now consider a finite momentum transfer with q = \q\ k = \k\ ^ I^fI- Due to the rotational 
invariance of the Fermi surface we can choose q to lie along the z axis, see Fig. 14, such that 

^c{k + q) vf{(I),0) ■ [k + q — kF{(t>,9)] = VF{k — kF) + vf{(I>, 9) ■ q = VF^k — kF) + vfQCOS^O) 

= ^dk) + VFqcos{9) (B20) 


which allows us to write the above sum of momenta as 

-I p^ItZ n\ poo 

^f[^v{k)dc{k +q)] = I d(j) dcos0 dk ■ k"^ f[-^cik)dc{k) + VFqcos{0)] 

^ Jo J-i Jo 

= / ^ J ^^^P^0f[-t^ + VFqcos{9)] 

(j(y r°° 

= + (B21) 

where B is the energy distance from the lower band edge to the Fermi energy in the conduction band. We will now 
explicitly calculate the imaginary part of (B18) for the symmetric particle and hole bands in the normal conducting 
state and omit the real part since it is featureless and gives just a small unimportant correction to the dispersion of 
the photons. We find for zero fermionic temperature Tp = 0 the imaginary part 

/ Virq 7 pOO / \ 
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' ' ■UFg 


(B22) 


The real part will be neglected in the normal state as it is only weakly energy and momentum dependent and gives 
no important renormalization of the photon excitation energies w* « Wq. 


2. Superconductor 


In the superconducting state we can use the standard parametrization of the normal and anomalous Greens functions 
{Tf = 0) 


/I-RM '^A,k,v/c 

^o,k,v/cy‘^) - ~ ^ TTbTTd + 


P A,k,v / c 


(B23) 


u! — Ey/c{k) ±i0 uj — Ey/c{k) ± iO 
K.k,v/ci‘^) = [l-2nF(w)] (A^,k,v/ci^) - ^ik,v/ci^)^ = -2TTi(^aA,k,v/c5[uJ - E^/^{k)] + PA,k,v/J[i^ + Ey/c{k)]^ 


where we defined the superconducting dispersion Ey/c{k) = \J^v/c{kY + |^d/c|^ and the coherence factors for the 
different Green’s functions are given by 
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An explicit calculation of the normal and anomalous self-energy show that 
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and the other two self-energies are related via the symmetries (B16). 


a. Anomalous self-energy 


Let us now consider first the anomalous self-energy that can only occur in the presence of superconducting 

quasi-particles. We can use the linearization (B21) for the symmetric bands to write the retarded self-energy as 
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where we assumed |Ac| = |A„| = A, defined w = w/A and approximated the DOS to lie near the Fermi surface since 
the dominant part of the integral will come from there. The real part of the integral is obviously symmetric and the 
imaginary part anti-symmetric in frequency 
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c\,R I 


imn-’"(o;) = -imnri“(-o;) 

such that we will continue to focus on positive frequencies w > 0. Let us first calculate the real part 
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for w > 0. We see that the integrand will have a singularity if 
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and we therefore expect that at w = 2 some important feature can occur. Let us now consider this case w « 2A, 
where it is convenient to separate the singular region Xmin ~ —a/2 of the integrand 
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(B30) 


The contribution gq{uj) will just give a small unimportant term as I^q] ^ A is small and therefore our dispersion 
is just changed slightly. The important contributions will come from the divergences at w = ±2A, which are captured 
well by the expression (B30), such that we approximate 
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2(^) ' 12-MI ' 

The imaginary part can be calculated using the same methods as just explained for the real part and is given by 
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Obviously, there is a gap of 2A for the imaginary part of the particle-hole bubble Imn'j\’^(a;) ~ d{uj — 2A). Let us 
again focus on the region w > 2A, where the dominant contributions come again from x « —a/2 such that 
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where in the end we approximated the integral for small a due to the integrand restrictions. Far away from the gap 
region w ^ 2 A it is easy to show that the imaginary part Imn“*|^^(cLi >> 2 A) ~ 1 /w^ falls down rapidly and is 
therefore not important, such that we can approximate Imn‘^*-|’^(a;) for all lu with (B33) as stated in the main text. 


b. Normal self-energy 


The normal self-energy around a;_ = uj — eVo = ±2A can easily be shown to behave exactly like the 

anomalous self-energy around uj = ±2A. Far away from this region |a;_| ^ 2A it just behaves like in the 

normal conductor (B22). In between, we interpolate between these two by replacing the self-energy as the maximum 
of the two limits 
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Appendix C: Derivation of the dressed propagators 


Let us first consider the Wigner transformation of a convolution of two bosonic matrices C = A o B with time 
dependence 
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where as usual we set (j){2T) = 2eVoT — (fc (fjv 2(j)g, defined T = and t = t — t'. The Wigner-transform of a 
convolution C = Ao B can be expressed conveniently by the Wigner transforms of A and B as 
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Using the identities 
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The Dyson equation in Keldysh space is formally given by the convolution 


(C5) 


T®1 _L Trbath 

advanced and Keldysh propagators of the photon 


where tig = fl®' + is the full self-energy of the system. The three coupled equations for the full retarded, 
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can now be solved with the ansatz that 
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have exactly the structure of A we can use (C4) to calculate the Wigner transform of the convolutions in (C6) exactly, 
which leads to the propagators defined in (54). 
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